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Abstract 

We perform a fit to precise electroweak data to determine the Higgs and top 
masses. Penalty functions taking into account their production limits are included. 
We find run = 65^4 45 GeV and mt = 122^20 GeV. However whereas the top x 2 
distribution behaves properly near the minimum, the Higgs x 2 distribution does not, 
indicating a statistical fluctuation or new physics. In fact no significative bound on 
the Higgs mass can be given at present. However, if the LEP accuracy is improved 
and the top is discovered in the preferred range of top masses, a meaningful bound 
on the Higgs mass could be obtained within the standard model framework. 
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Introduction 



The standard model of electroweak interactions [I] describes impressively well all experi- 
mental data 0. This description depends on two unknown parameters, the top mass, m t , 
and the Higgs mass, mg. The present experimental accuracy of the electroweak measure- 
ments justifies to look for indirect bounds on these masses. In this paper we present the 
X 2 distributions resulting from the most precise and recent electroweak data as functions 
of m t and m H . Similar analyses can be found in the literature £|, |[. We pay particular 
attention to the Higgs penalty function describing the Higgs production limit. It plays 
an important role due to the low sensitivity of the data on m# and to the fact that the 
minimum of the x 2 i if only indirect data are used, is (at present) well below the direct 
production limit. Actually, the present sensitivity of the data on itlh will prove to be too 
low to obtain stringent, statistically significant, limits. However the sensitivity could be 
greatly improved in the near future, provided that LEP errors get reduced and the top is 
discovered. 



Data 

The experimental values used in the fits are gathered in Table 1. The detailed way LEP 
excludes Higgs production can be found in Ref. ||. For our discussion it is enough to 
realize that the final result of the search is that no event with a given signature is found, 
whereas some events should be present if the Higgs boson would have been produced. Of 
course, the number of expected events decreases with increasing m#. Therefore, assuming 
a Poisson distribution for the number of expected events, values below m# can be excluded 
(for instance) at the 95% ci if the number of expected events for m# is larger than or 
equal to 3. The number of expected events at LEP, N, as a function of rriH is shown in 
Fig. la (this plot is obtained from the compilation in Ref. J7|). The confidence level, ci, 
corresponding to seeing no events when N events are expected is equal to 1 — e~ N ; and 
this confidence level is transformed into the corresponding number of standard deviations, 
x, for a single sided N(0, 1) distribution using the relation 



c£ + 1 



/X „ 
e~* dt. 
-oo 



In this way x can be written as a function of m#. This function is plotted in Fig. lb, and 
the best linear fit, 

x = 15.81 - 0.24 m H , 

is superimposed on it. Thus, the corresponding contribution to x 2 is given by the penalty 
function 

f [ 65.0-m H (GeV) \ 2 
Ax 2 H =\ { Tl ) Jm H <65GeV, 

{ if m H > 65 GeV. 

In the case of the top quark the detailed way TEVATRON excludes top production 
can be found in Ref. pi. In this case, however, the minimum of the x 2 ■, is we U above the 
direct production limit, and then the corresponding penalty function plays no important 
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role in the fit. The strategy for the top search is also different. What is studied in 
the top case is the presence of an excess of events in some distributions (basically in 
the transverse mass distribution) due to top production. These distributions are fitted 
with a linear superposition of the standard model predictions with no top quark plus 
the top production. The result of this fit for each top mass is a probability density 
distribution L(p) for the variable p defined as the ratio of the number of fitted to the 
number of expected (a t (m t )) top events for such a mass. No top events corresponds to 
p = 0. To exclude a top mass (for instance) with a 95% cl we associate to each mass 
m t the value of po{mt) which defines for the probability density corresponding to this 
mass an area, J po(mt) L(p)dp, equal to the 95% of the total area or probability, Jq°° L(p)dp. 
Then po(m t ) = 1 gives the m t limit below which the top mass is excluded. To make 
this computation for different confidence levels we have assumed that p (mt)at(m t ) is 
independent of m t in the region of interest. (Note that this is equivalent to assume 
that the shape of the top distribution is independent of the top mass in this region.) 
To translate the ci to which a given top mass is excluded into a number of standard 
deviations we proceed as in the case of the Higgs mass. The corresponding fit gives the 
approximate penalty function 

Ax?= { 9 J ifm t <108GeV, 

[ if m t > 108 GeV. 

The x 2 a l so includes neutral current data: neutrino-quark data, vq ||, the latest 
neutrino-electron results, z/„e [ffD], and the parity violation in atoms data, eH JTTI, O; the 



W mass limits: M w |L3| and [|14]]; the LEP results presented at the Joint International 



Lepton-Photon Symposium and Europhysics Conference on High-Energy Physics in Au- 
gust 1991 (with the correlation matrix given in Ref. [[HJ); and the strong coupling 



constant, a s WA - We take for a s the latest (most precise) ALEPH measurement using 
hadronic Z decays. The definite a s value is important to fix the x 2 value at the minimum 
but modifies little the relative x 2 distribution and does not affect the conclusions. 



Theory 

We use the on-shell scheme of Ref. [|18[] with the electromagnetic coupling constant, a, and 
the Fermi constant, G^, as input parameters. For the calculation of the observables in 
the Z physics sector two independent electroweak libraries have been used: the one from 
G. Burgers and W. Hollik and the one from the Dubna-Zeuthen group | 20| . In spite 



of the fact that they use a completely different calculational scheme, their results are in 
very good agreement |[21||. We have upgraded the first one to incorporate: 



The missing relevant pieces of two- loop corrections p2| ; 
The dominant QCD corrections to m t dependent terms 



The updated calculation of the QCD corrections to the decay Z — > bb [24|; and 
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0(al) corrections to the decay widths 



25 1 



Similar upgrades have been implemented in the Dubna-Zeuthen library by their authors. 
A comparison of the present versions of both programs shows an even better agreement 
than the one previously quoted. While the numbers and plots shown in this study have 
been obtained with the first library, we have explicitly checked that the changes are 
negligible if the Dubna-Zeuthen library is used instead. In our fits we use the explicit 
expressions for all the observables, in particular for the neutral current ones |26| . We have 
checked that reducing the full experimental information for the latter to the equivalent 
electroweak mixing angle value, our fits and conclusions change little. 



Interpretation 

Our results are summarized in Figs. 2-8. Fig. 2a shows the contribution to the x 2 distribu- 
tion of neutral current and pp collider data; Fig. 3a that of LEP data and a s ; Fig. 4a their 
sum; and Fig. 5a that of all the data including the penalty functions. Thus Figs.2a-4a 
give partial contributions to Fig. 5a. At each point in the m t — rriH plane we minimize 
with respect to M z and a s (G M and a being fixed). Thus, the total number of degrees of 
freedom is 20 — 2 = 18. (The two degrees of freedom corresponding to the penalty func- 
tions are taken care by the two parameters, m t and tuh, with respect to which we do not 
minimize.) Figs.2b-5b show the same fits but replacing the central experimental values 
of the different observables by their predicted (standard model) values near the minimum 
of the x 2 distribution. In particular we have taken m# = 70 GeV, m t = 130 GeV. The 
corresponding standard model predictions for the observables used in the fits are given in 
Table 2. Comparing Figs.2a-4a and Fig.5a we see that the Higgs penalty function moves 
the minimum from m# = 10 to 70 GeV. However, the main observation to be elaborated 
below when comparing Fig. 3a and Fig.3b is the apparent change in the LEP distribution, 
and then in the global one. This tells us that some experimental central values manifest 
a statistical fluctuation or indicate new physics. This comparison also proves the actual 
lack of sensitivity of the data on ran. To understand this change it is illuminating to 
rewrite (in obvious notation) 

, _ (XT - X*(m H , m u ...)f _ (XT ~ X? + e,f 
X{m H ,m t ,...) — 2 ~ 2-i 



at " of 



E 



(Xr - X»f - 2(Xr - X^ + e\ 
of 



where X® are the standard model predictions at the x 2 minimum and 6j are the differences 

between the predicted values at a given point and the values at the minimum. The 

last expression shows that a statistical fluctuation or a relatively large change in the 

central value of any observable contributes to the x 2 with a relatively large constant term 

(XT - X°? j .,1 y , , , . -2(AT P -X°)e t 

2 an d with a linear term of pronounced slope ^ • However, 

® i @i 

e 2 

what gives the sensitivity of the observable is — ^. Comparing Figs.2a-5a and Figs.2b-5b 

erf 
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we see that indirect LEP data are very insensitive to van but some observables have a 
large discrepancy between their measured central values and the predicted ones. The 
global fit to all data (Fig. 5a) gives 

m H = 65lf 5 GeV, 
m t = 122+^5 GeV. 

Although the fits in Figs.2a-5a give the actual experimental information, those in Figs. 2b- 
5b give the actual sensitivity of the data on mn and m t . Hence present indirect elec- 
troweak data are not sensitive to m#. Actually, the present LEP discrepancy can be 
traced back to the value of the ratio Re and the bottom forward-backward asymmetry 
A b p B . This is to say that the data show a higher sensitivity than expected on m H due to 
the pronounced slope terms (introduced above) resulting from Re and Ap B . If this dis- 
crepancy is a statistical fluctuation and it disappears when more data are available, the 
subsequent fit will approach Figs.2b-5b. However, it can also happen that the central val- 
ues of other observables fluctuate, in which case our discussion would apply to them. At 
any rate this will be always the case (if new physics is discarded) when the preferred m# 
value from indirect LEP data lies below the production limit. In Figs.6a,b we show the \ 2 
distribution as a function of m# (minimization with respect to m< is understood) for the 
global fits in Figs.5a,b. Plotting the same distribution for the different observables it can 
be explicitly proven that only Re and A b p B show the pronounced slope behaviour for the 
real data, as can be seen in Fig. 6a, whereas the other observables show a small variation, 
which we do not plot. Fig. 6b confirms that the sensitivity of Re and Ap B on m# is very 
low, as it is the sensitivity of the other observables and of the full set. In fact we plot A% 2 
in Fig. 6a, subtracting from each x 2 contribution its value at the global minimum: 10.79, 
0.88, 2.73 for the global, Re, A b p B contributions, respectively. Using higher values for a s 
would arrange the Re deviation from the standard model prediction. This would translate 
into a modification of the mt value at the minimum in order to conserve the agreement 
between the measured and the predicted I 'z values. However A b p B would still constrain 
the fit and keep unchanged the value of m# at the minimum and the conclusions. Lower 
values of a s make worse the disagreement between the measured and the predicted Re 
values. It is worth to note that although the top mass limit is significative, it is strongly 
correlated to the Higgs mass as can be seen in Fig.5. For instance if the Higgs mass were 
known to be 300 ± 30 GeV, the present data would imply m t = 139^25 GeV (see Fig. 5a). 



Conclusions 

As stressed above, today's data show a fictitious m# sensitivity due to a statistical fluctu- 
ation or to new physics manifested by the experimental values of Re and A b p B (Figs.2a-6a). 
As a matter of fact, the actual sensitivity of the considered observables on m# is very 
low (Figs.2b-6b). The Higgs penalty function plays an important role in the fits, push- 
ing up the Higgs mass value at the minimum. What is more interesting is the expected 
improvement of the m# sensitivity in the near future | S7| . Figs. 2-5 show the correlation 



between wh and m t which should translate, once the top quark is discovered and its mass 
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determined, into a definite bound on ran- In fact, if present LEP accuracy is improved as 
assumed in Table 2, and the top mass is in the preferred range given above, a meaningful 
bound on the Higgs mass (within the standard model framework) could be obtained. This 
fact is illustrated in Figs. 7 and 8. In Fig. 7 we plot the ^-distribution using the standard 
model predictions as central values and the set of LEP improved errors quoted in Table 
2. We minimize at each point with respect to Mz and a s . Fig. 7 shows that fixing m t a 
relatively strong upper bound on itlh can be deduced. This is explicitly shown in Fig. 8, 
where we fix the top mass, m t = 130 ± 1 GeV. In this case m# < 315 GeV at 95% ci. 
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Table captions 

Table 1 Experimental values of the observables used in the fits. 

Table 2 Standard model predictions for the observables used in the fits and for ran = 
70 GeV and m t = 130 GeV. The improved error column gives a guess of future 
LEP errors. 



Figure captions 

Fig.lafb] Number of expected events [standard deviations] at LEP as a function of m#. 

Fig.2a[b] Level contours for the contribution to the ^-distribution of neutral current 
and pp collider data in Table 1 [standard model predictions in Table 2 with the 
present experimental errors in Table 1] . Minimization with respect to Mz and a s is 
understood. 

Fig.3a[b] The same as in Fig.2a[b] but for LEP observables including a s . 

Fig.4a[b] The same as in Fig.2a[b] but for all observables in Table 1. 

Fig.5a[b] The same as in Fig.2a[b] but for all observables in Table 1 plus the top and 
Higgs penalty functions given in the text. 

Fig.6a[b] ^-distributions as functions of mu for the global fit in Fig.5a[b] (solid curve). 
The Re (dashed curve) and the Ap B (dotted curve) contributions are also shown. 
Minimization with respect to Mz, a s and m t is understood. 

Fig. 7 Level contours for the total ^-distribution assuming the standard model predic- 
tions as central values and the set of improved errors given in Table 2. Minimization 
with respect to M z and a s is understood. 

Fig. 8 ^-distribution as a function of m H for the global fit in Fig. 7 and m t = 130±1 GeV. 
Minimization with respect to M z , a s and m t is understood. 
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Quantity 


Experimental Value 


Correlation Matrix 


vq 


9l 
9r 
6 L 
Or 


0.2977 ±0.0042 
0.0317 ± 0.0034 
2.50 ±0.03 
4.59+S-tl 

— U.z ( 


1. 


1. 


1. 

1. 




9a 
9 e v 


-0.503 ±0.018 
-0.025 ±0.020 


1. 


-0.05 
1. 




eH 


Ciu + C\d 
Ciu — C\d 

Clu — Cid 


0.144 ±0.007 
-0.60 ± 0.09 
-0.05 ±0.11 


1. 


1. 


1. 


pp 


M w 
M w /M z 


79.91 ± 0.39 GeV 
0.8813 ±0.0041 






M z 

< 
Re 


91.175 ±0.021 GeV 
2487 ± 10 MeV 
41.36 ± 0.23 nb 
20.92 ±0.11 


1. 


0.09 
1. 


0.01 0.00 
-0.25 -0.07 
1. 0.18 
1. 


LEP 


From Afb- 

( 9 i(M Z )) 2 

\9a J 


0.0048 ±0.0012 


1. 








From A pol (r): 

9 f(Mz) 
9a 












0.072 ±0.017 




1. 






From 66-asymmetry: 
A%(M Z ) 


0.132 ±0.022 






1. 




From gg-asymmetry: 
sm 2 e w (M z ) 


0.2303 ± 0.0035 






1. 




a s 


0.125 ±0.005 





Table 1 
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Quantity 


Standard model prediction 


Improved errors 


uq 


9l 
9k 

0L 
0R 


0.3012 
0.0302 
2.46 
5.18 


0.0042 
0.0034 
0.03 
0.44 




9a 
9v 


-0.505 
-0.048 


0.018 
0.020 


eH 


Ciu — C\d 

C2u — C2d 


0.146 
-0.54 
-0.09 


0.007 
0.09 
0.11 


pp 


M w 
M w /M z 


80.14 GeV 
0.8790 


0.06 GeV 
0.0041 




M z 

T z 

o° h 
Re 


91.175 GeV 
2489 MeV 
41.41 nb 
20.81 


0.005 GeV 
4 MeV 
0.08 nb 
0.02 


LEP 


From Afb- 








( 9 i(M Z )) 2 


0.0048 


0.0006 




From A pol (r): 

p 








air 

9 -i(M z ) 
9 a 


0.069 


0.009 




From 66-asymmetry: 
A%(M Z ) 


0.097 


0.010 




From gg-asymmetry: 
sin 2 6 W {M Z ) 


0.2320 


0.0035 




a s 


0.125 


0.005 



Table 2 
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